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Abstract. Using the norm as a measure of mixing, we prove that 2d 

Euler flows on the torus mix passive scalars at most exponentially. The mixing 
rate is bounded linearly by the BMO norm of the vorticity (and thus by its 
L°° norm). We also give an analogous bound on the growth rate of scalar 
gradients. 



1. Introduction 

Let v(x,t) be the velocity field of an incompressible inviscid fluid flow, governed 
by the Euler equation 

(1.1) dtv + v-Vv + Vp = with V-v = 0. 

For concreteness, we take {x,y) ^ x E Ai := [0, 1]^ with periodic boundary con- 
ditions. With no loss of generality, we assume that v = (u, v) has zero integral 
over the domain M. For our purpose here we assume that the initial data v(-,0) 
is smooth, say, C^. The evolution of a passive scalar 6{x,t) in this flow is then 
governed by 

(1.2) dte + v-ve = o. 

As with V, we assume that the integral of 9 over M vanishes. 

Assuming the boundedness of v for all t > 0, one can show from ()1.2p that \9\i[^p 
is constant in time for all p £ [l,c»]- It is clear that there are many solutions of 
(|l.ip . e.g., shear flows and (nearly) stable ones, for which the scalar 6 will undergo 
little time evolution. If the solution of (|l.ip is chaotic (as is often the case in 
many interesting situations), however, one expects that 9 will become increasingly 
"mixed" in time; mathematically, one expects that 9{-,t) will converge weakly to 
as i — > GO in some suitable space. Of significant mathematical and physical interest 
is the rate of this mixing or convergence. 

In [9], the authors introduced a coarse-grained "mix norm" which they showed 
to be equivalent to the H~^/^ norm, defined formally as 

(1.3) \w\%.:=j:^^,\k\^^\wk\^ 

where Wk are the Fourier coefficients of w. (This is a seminorm in general, but is 
a proper norm when w has zero integral over A^, which will always be the case in 
this work.) They argued that this H~^^^ norm is a good measure of mixing as it is 
weighted towards larger-scale features at the expense of fine details. 
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(1.7) -\Vip\l2 > -ci \Vv\L^\Vip\l.. 



Subsequently, [8] proposed the use of the more convenient H^^ norm to measure 
the mixing rate, having proved that for any s < and for any 0[-,t) G L'^ with 
zero integral, the weak convergence d{-,t) in as i — > cxd is equivalent to 
|0(-,t)|^, ^0. 

In this work, we follow their use of the H ^ norm as our measure of mixing. 
Putting if and multiplying ()1.2p by —tf in L^, we have 

(1-4) ^^l^'^li^ - {v-V/\ip, ^)^. - 0. 

Assuming sufficient smoothness, we integrate the nonlinear term by parts 

Using the standard estimate 

(1.6) \{{d,v)-Vip,d,f)^,\ <c \Vv\L^\Vip\l2, 

we have 

d_ 

dt' 

Without assuming that v is a solution of (|l.ip . the possibility of perfect mixing 
in a finite time T was not ruled out in [5], where lower bounds for perfect mixing 
time were given for physically relevant cases; here perfect mixing is defined by 
\e{-,T)\j^_i = \Vip{-,T)\^2 = 0. They did however point out that if \yv{-,t)\^^ is 
uniformly bounded, this will give a bound on the mixing rate. 

When V is the solution of (jl.ip . however, one can use the argument of [3 [S] to 
rule out perfect mixing as long as the solution remains regular. Unlike the results 
in the following sections, this argument works both in 2d and 3d. We fix s > 3 
and assume that v{-,0) G H^. By Sobolev embedding, this implies Wv{-,0) G L°°. 
Using the differential inequality 

(1.8) ^^\v\h <(^2\^v\L^\'^\h, 

we have 

(1.9) \v{;t)\l. < \v{;0)\l.e^v(^C,j\\7v{;t')\^^ dt') . 

Now (II. 7p can also be integrated to give 

(1.10) |V^(-,i)|i. > |V^(-,0)|i.exp(-Ci / \Vv{;t')\L^dt'). 



In the 3d case T], one can bound |Vt>|^oo essentially by jcurlvj^^oo log \v\^3 to show 
that the integral in (jl.9p is finite for alH > such that curlt; G L^{[0,t]; L°°). In 
other words, no perfect mixing is possible before the blow-up time if the latter is 
finite. In the 2d case, it is shown in ^ that v{-,t) G H'' for all t > 0, implying 
that the integral in (jl.lOp is finite and | V(y3(-, t)|^2 > for all t > 0. Therefore, 
no perfect mixing is possible in finite time when is a solution of the 2d Euler 
equation (|l.ip . 
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2. A Uniform Bound on the Mixing Rate 

The inequality (jl.lOp does not tell us much about the mixing rate since we do 
not have a bound uniform in t for \Vv\]^aa- As noted above, [T] bounded |Vv|^oo 
essentially by |curli;|^oo log |7;|j:^3, which is not bounded uniformly in t — in fact, 
one expects \v{-,t)\jj3 to grow without bound for "generic" flows. 

In this work, we take advantage of the special structure of the Jacobian to bound 
(|1.6|) by |Vv|bjjo instead of iVvj^^cx, and use a time-uniform bound on |Vv|j3j,,() 
to obtain a bound on the mixing rate. To this end, we introduce the vorticity 
u :— dxV — dyU, whose integral over Ai vanishes by Stokes' theorem, and the 
streamfunction ip := A~"'^a; (defined uniquely by requiring that its integral over Ai 
vanishes), in terms of which v = V-^ip = {—dyip,dxip)- Taking the curl of (jl.ip . we 
have 

(2.1) dtLu + diij,io) ^0 

where the Jacobian d{i/j, w) := V^-^'^"^ = dx'ip dyU — dx^-o dyi/j. Using (|1.5p . we can 
write (jl.4|) as 

(2.2) ^^|v^ii.+E, (9(9,V;,^),a,^)^. -0. 

For our purpose here, the following definition suffices; we refer the reader to [12J 
for more details. Let A/" C M^. For w G Ll^^{M) and any ball B C TV, let denote 
the average of it; in B, 

(2.3) Wb ■■= J^w{x) dx. 

We say that w € BMO (TV), the space of functions of bounded mean oscillations in 
A", if 

(2-4) ^Ibmo := sup [ \w{x) ~ Wb \ dx < oo. 

We note that | • |b^,jq as defined here is a seminorm (constants have BMO-norm 
zero), but since we only deal with functions of zero average, (j2.4l) defines a proper 
norm for all relevant quantities. We also note that (12.41) implies 

(2.5) h|BMO<klL~- 

Our main result is the following: 

Theorem. Let the passive scalar 9 be the solution of (jl.2p and v that of the Euler 
equation (jl.ip with initial vorticity a;(-,0) G L°°[A4). Then 6 satisfies 

(2.6) \6{;t)\l-, > \e{-Ml-^ "^""pI"-^/* H-^t')Uo dt') 

for some constant X depending only on the domain M. Moreover, the mixing rate 
is bounded by the sup-norm of the initial vorticity as 

(2.7) \e{;t)\l^, >|0(.,O)|?,-,exp(-a|c.(.,O)L.o). 

We note that our result does not address the more difficult issue of the existence 
of exponentially mixing solutions of (jl.ip . or whether the bound (|2.7I) is attained 
even qualitatively. See [1] and references therein for further discussion. 
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Proof. Parts of this proof were inspired by As usual, c denotes a generic positive 
constant whose value may differ each time the symbol appears. Assuming for now 
the following estimate for the Jacobian, 

(2.8) |5(C,^)L. <c|VC|bmo|V¥'L2, 

we take = dxip and ( — dy^p in turn in (|2.2|) to get [cf. (|1.7p ] 

(2.9) ^iV^li. >-c|Vt;|3,,jV^|i.. 

This implies the analogue of (jl.lOp with BMO in place of L°° . Unlike the L°° case 
where |w|^oo does not bound |Vt>|^oo, here we have 

(2.10) |Vi;|,„„<c|w|3,„, 

whose proof will follow shortly. Using ()2.10p in ()2.9p and integrating gives us p.6p . 
Using ()2.5p and the fact that \uj{-,t)\]^^ ~ |a;(-,0)|^oo gives us <\2.7^ . 

To prove (|2.10p . we start with the identity (cf. 11, p. 59]) 

(2.11) Vv = VV^A-^cj = {R,,Ry){Ry, -i?,) u 

where '■— dx{—A)^^/^ and analogously for Ry. It is shown in \TT, p. 138] that 
the Riesz transforms {Rx,Ry) are bounded operators in the Hardy space ■H^(K^), 
which is dual to BM0(R2). Boundedness of R := {Rx.Ry) in BM0(R2) then 
follows by duality. From the definition (|2.4p . it is clear that the norm in BMO(A^) 
is identical to that for periodic functions in BMO(M^). 

We now prove the Jacobian estimate (|2.8p . For bounded domains in with 
Dirichlet boundary conditions, although not stated explicitly, the proof is essentially 
contained in [6]. Its extension to periodic boundary conditions is essentially done 
in [5]. We reproduce these proofs here (with minor modifications) for convenience. 

Let M.I := Ux^mBi{x) where Bi{x) C R^ is the ball of unit radius centred at 
X. We first prove ^ for C G H^iMi) D {VC e BMO(Xi)} and if G H^iMi). 
Denoting the Fourier transform by a hat, we use the following result from O p. 154]. 

Let ct(-, •) G C°°(R2 X R2 - {0}) satisfies 

(2.12) |9|'9,^a(e,r,)| <a;3(|C| + M)-H-l^l V^, r; G R^ - {0} 
for every multi-indices a and /?, and 

(2.13) a(0,?7) = 0. 
Then the bilinear operator 

(2.14) a{f,g){x):^ ( e'-<^+^K{^,v) fiO gin) dvd^ 

is bounded in L'^{A4i) as 

(2.15) \<j{f,9)\L^ <c\fU,\gy. 

For the Jacobian 9(C, f), we have i7(^, r]) = — • r], so in order to satisfy (|2.12l) - 
(|2.13p we split d{C, ^p) as follows. Let p G C°° be monotone increasing with p{t) = 
for i < 1 and p{t) = 1 for t > 2, and let d{C, f) = J> + J< where 

J>{x) = -^ f e'"-(«+'')p(|C|)(e^-r,)C(O^WdrydC 

(2.16) 

J^i^) = -7r / e'-(«+'')[i-p(iei)](e^-r?)c(e)^(^)d,7de 
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Taking — p(|^|), we have from ()2.15p 

(2.17) |J>L. <c|VC|bmo|V<^|l2. 
And taking 

(2.18) ^(?,.)-^Y^e\ 
(|2.15p gives us 

(2.19) |J<L. <c|VC|bmo|V^|l.. 

Having proved (j2.8p in A4i with Dirichlet boundary conditions, we now extend 
it to the periodic case, taking ( G i?por(R^) H {VC £ BM0(M2)} and tp e H^^^X^"^). 
Let g e C°°(R2; [0, 1]) with g = 1 in and g = in - Xi. By the above, we 
have 

(2.20) |a(^.C,ML.(^,) <C•(^^,>^l)|V(^.C)lBMo(>.,)l^(^'^)li^(Ml)■ 

By the periodicity of C and smoothness of g, we have 

(2.21) |V(eC)lBMo(>i,) < C{M,Mi,g) |VC| 
and similarly 

(2.22) |V(M|i.(^^) < C{M,Mi, g) 

One then takes the infimum over g in (|2.2ip - (|2.22l) to remove the dependence of 
the constants on g. The conclusion follows from these and 

(2.23) \d{C,v)\L2(M) = I^(£'C,£«P)L2(^) < |5((?C,e"/')L2(^^)- 

□ 

3. Further Remarks 

One can obtain an upper bound on the scalar gradient using the same method: 
Multiplying by -A^ and estimating, we find [cf. ([T3t and (^3)) ] 

(3.1) Z^'^^''^ = ^'^^^ {{d,v)-ve,d,e)^. 

<C(A^)|Vt;Uo|V0|i.. 

As before, since |Vij|b^,,o < c|w|^oo and using the fact that the latter quantity is 
time-invariant, we have [cf. (|2.7p ] 

(3.2) \Ve{;t)\l, < \Vd{;0)\l.e^p{tX{M)\uji;0)\^^). 
Since the vorticity equation (12.11) . or equivalently, 

(3.3) dtLJ + v-Vuj = 0, 

is of the same form as (|1.2p . these bounds evidently apply equally well to j Va;(-, i) |^2 • 
For this active scalar, one expects Vw to grow without bound; see, e.g., [lOj for 
estimates in Holder spaces. 



Acknowledgment. The author thanks Jacques Vanneste and Xiaoming Wang for 
insightful discussions leading to this work. 



D. WIROSOETISNO 



References 



[1] J. T. Beale, T. Kato, and A. J. Majda, Remarks on the breakdown of smooth solutions 
for the 3-D Euler equations, Comm. Math. Phys., 94 (1984), pp. 61-66. 

[2] R. R. COIFMAN AND Y. Meyer, Au dela des operateurs pseudo-differentiels, Asterisque, 57 
(1972), pp. 1-185. 

[3] S. Gottlieb, F. Tone, C. Wang, X. Wang, and D. Wirosoetisno, Long time stability of 

a classical efficient scheme for two dimensional Navier-Stokes equations, I arXiv: 1105 .4349 
[4] P. H. Haynes and J. Vanneste, What controls the decay of passive scalars in smooth flows?, 

Phys. Fluids, 17 (2005), p. 097103. 
[5] T. Kato, Remarks on the Euler and Naviei — Stokes equations in R^, Proc. Sympos. Pure 

Math., 45;2 (1986), pp. 1-7. 
[6] N. Kim, Large friction limit and the inviscid limit of 2d Navier-Stokes equations under 

Navier friction condition, SIAM J. Math. Anal., 41 (2009), pp. 1653-1663. 
[7] H. KOZONO and Y. TaniuCHI, Bilinear estimates in BMO and the Navier-Stokes equations. 

Math. Z., 235 (2000), pp. 173-194. 
[8] Z. Lin, J.-L. Thiffeault, and C. R. Doering, Optimal stirring strategies for passive scalar 

mixing, J. Fluid Mech., 675 (2011), pp. 465-476. 
[9] G. Mathew, I. Mezic, and L. Petzold, A multiscale measure for mixing, Physica, D 211 

(2005), pp. 23-46. 

[10] A. MORGULIS, A. Shnirelman, and v. Yudovich, Loss of smoothness and inherent insta- 
bility of 2d inviscid fluid flows. Comm. P.D.E., 33 (2008), pp. 943-968. 

[11] E. M. Stein, Singular integrals and differentiability properties of functions, Princeton Univ. 
Press, 1970. 

[12] , Harmonic analysis, Princeton Univ. Press, 1993. 



E-mail address: djoko.wirosoetisno(3durh5mi. ac.uk 
URL: http : //www .maths . dur . ac . uk/"dinaOdw 

Mathematical Sciences, Durham University, United Kingdom 



